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ABSTRACT 

The  derivation  of  the  stress-integral  form  of  the  Higashitani-Pritchard-Baird-Lodge 
relation  cannot  be  applied  to  plane  Poiseuille  flow.  It  ignores  the  contribution  of  the 
streamwise  pressure  gradient,  and  doing  so  can  lead  to  contradictions.  Furthermore,  the 
variable  change  which  leads  to  the  stress-integral  form  of  the  HPBL  relation  is  not  valid  in 
a  second-order  fluid  approximation  for  slow  flows.  Correcting  these  deficiencies  leads  to  a 
modified  relation  which  appears  intractable.  Nevertheless,  the  Tanner-Pipkin  result  that 
Pe  =  jVj/4  in  a  second  order  fluid  is  valid  when  properly  interpreted.  > 
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ON  HOLE-PRESSURES  IN  PLANE  POISEUILLE  FLOW  OVER  TRANSVERSE  SLOTS 

D.  S.  Malkus  and  Minwu  Yao 


1.  Introduction 

The  symbols  used  in  this  paper  are  defined  in  Appendix  I  at  the  end  of  main  text. 
The  problem  in  question  is  plane,  creeping  flow  over  a  transverse  slot,  which  is  illustrated 
in  Figure  1  and  described  in  refs.  1-4.  Ref.  1  describes  a  proposed  relation  between 
pressure  differences,  Pe.  and  the  first  normal-stress  difference,  JVj.  It  is  clear  that  ref.  1 
is  primarily  concerned  with  Couette  flow  (plate-driven)  and  not  Poiseuille  flow  (pressure- 
gradient  driven).  In  refs.  2-8  the  possibility  of  applying  the  relation  to  Poiseuille  flow 
has  been  explored  experimentally  and  numerically.  The  final  relation  is  attributable  to 
Higashitani.  Pritchard,  Baird,  and  Lodge  [2,3] ,  and  has  two  forms  for  creeping  flow: 


_  rflnP, 

-  ’  1  U'  _  ,  r  f 

alnou. 


where  ow  and  Niw  are  wall  values  of  the  primary  viscometric  functions,  and  Pe  is  the 
pressure  difference  from  top  to  bottom  of  the  slot.  Eq.  (lb)  is  easily  derived  from  eq. 
(la),  and  the  purpose  of  this  paper  is  to  show  that  the  derivation  of  eq.  (la)  is,  at  best, 
inapplicable  to  plane  Poiseuille  flow.  It  is  with  some  irony,  then,  that  the  authors  can 
state  that  numerical  results  seem  to  corroborate  eqs.  (la)  and  (lb)  to  within  about  10%, 
at  worst,  for  a  Maxwell  fluid  and  to  validate  eq.  (lb)  to  u  working  rheological  accuracy” 
for  some  more  realistic  fluids.  These  results  will  be  discussed  thoroughly  elsewhere  i 7.8 ; . 
It  seems  to  be  an  unavoidable  conclusion,  however,  that  the  numerical  and  experimental 
results  owe  their  pleasing  concordance  to  some  other  cause  than  the  validity  of  the  stress- 
integral  form  of  the  HPBL  relation  eq.  (la),  at  least  as  it  is  currently  understood. 

This  paper  does  not  directly  concern  numerical  results,  but  it  does  often  rely  on 
evidence  which  can  be  observed  from  the  relatively  simple  problem  of  plane  Stokes  flow 
over  a  transverse  slot.  Since  the  computing  technology  at  the  time  of  this  writing  allows 
many  researchers  in  fluid  mechanics  to  obtain  and  analyze  solutions  to  Stokes  flow  problems 
with  relative  ease,  rather  than  tabulating  data  to  support  our  statements  about  Stokes  flow , 
we  invite  the  skeptical  reader  to  verify  them  directly.  We  believe  that  most  of  our  claims 
about  Stokes  flow  are  reasonable,  if  not  obvious,  and  should  not  provoke  controversy.  The 
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FIGURE  1 

Domain  and  dimensions  for  plane  flow  over  a  transverse  slot. 
Rheological,  Cartesian  co-ordinates,  x,,  have  origin  at  inter¬ 
section  of  slot  (£  and  centerline  of  channel  or  ‘die.’ 


reason  that  Stokes  flow  is  central  to  what  we  present  here  is  that  our  observations  are 
based  on,  but  not  confined  to,  the  second-order  fluid  range  of  non-Newtonian  behavior. 


2.  Preliminaries 


It  will  serve  the  purposes  of  completeness  and  facilitate  the  comparison  of  the  the¬ 
oretical  results  obtained  here  with  numerical  results  discussed  elsewhere  if  we  begin  by 
deriving  the  equation  for  plane  Poiseuille  flow,  taking  particular  care  to  fix  the  definition 
of  the  hydrostatic  pressure  function  in  the  same  manner  as  in  the  numerical  model  of  [4]. 
This  also  fixes  the  arbitrary  constant  in  p. 


2.1  Pressure  definition: 


In  plane  flow,  the  pressure  is  defined  so  that 


Ar,/2  Oy,  0 

o  12  —  A  i  '  2  0 
0  0  0 


'■>  1 


m 


mm 


where  A'j  =  on  -  a22.  In  a  shearing  flow.  A'i  is  the  usual  viscometric  function.  Eq.  (2) 
implies  that 

P=~{°  n— 022)/2  (3) 


In  addition,  we  specify  that 


pdV  =  0 


The  equations  of  motion  for  undisturbed  flow  (d  =  0), 

^n,i  ■+■  2,2  ~  0 


Ol\,l  +  ^22.2  =  0 


imply  that  for  plane  Poiseuille  flow 


^12, 2  -  Pa~  (W/2),,  -  k 
0 

OjA,\  -  (Ai/2),2  =  P> 2 

or,  imposing  a  boundary  condition  at  the  wall  and  symmetry  condition  at  mid-channel 

r  Q\i  —  kx 2 
-  uj(-,/i/2)  =  0 
.  u1>2(-,0)  =  0 

where  o12  =  r)^  and  ->(0)  =  0.  Integrating  eqs.  (6)  implies 

p  =  fcx,  +  /(x2) 
p  =  -AT,  /2  -  g(xi) 

thus 

frxi  -r  iVi(x2)/2  *  /(x2)  -  g(x i)  =  0 


From  eq.  (4) 


p  =  kx j  -  Afj/2  +  C 


1  fh/2 

-  /  ATi(x2)dx2 
J-h.l1 


and  for  a  Maxwell  fluid,  an  Oldroyd  B  fluid,  any  fluid  for  which  Arj  =  I'l'y2.  u \  -  constant, 
and  77  =  constant, 

C  =  A’lu,/6  (8) 


Eq.  (8)  holds  approximately  for  any  fluid  which  has  a  second-order  fluid  range,  when  the 
flow  is  sufficiently  slow.  For  such  fluids  and/or  flows  and  the  above,  we  may  deduce  the 
form  of  the  stress  tensor  in  undisturbed  plane  Poiseuille  flow: 

—  kx  i  —  A  i  —  ;\  j  u.  /  6  <ru.  0 

o  -  au  -kx\  -  Arju,/6  0  (9) 

0  0  —  kx  j  -f-  Ni /2  —  A^ity/6 

It  should  be  kept  in  mind  that  Niw  is  a  constant.  Observe  that  a  i2  =  aw  at  the  top 
wall,  as  illustrated  in  Figure  2.  For  the  pictured  flow  direction  ow  <  0.  A  more  crucial 
observation  is  that  0 22  is  a  a  linear  function  of  tj  alone  and  has  slope  -k;  this  will  soon 
become  important  to  our  arguments.  o\\  varies  across  channel  as  a  function  of  N\. 


FIGURE  2 

Undisturbed  plane  Poiseuille  flow,  showing  stress  tensor  at  the 
wall  where  0\2  =  +ow  (a  negative  value  for  the  pictured  flow). 
oX2  —  —Ow  at  the  opposite  wall,  and  the  stress  tensor  is  oth¬ 
erwise  the  same.  Can  be  thought  of  as  slot-flow  with  a  d  —  0 
slot. 


3.  Tanner-Pipkin  Results 

We  shall  primarily  be  concerned  with  slow  flows,  in  which  for  many  constitutive 
equations,  the  second-order  expansion  is  valid  9.10  .  In  such  cases  we  presume  that  by 
solving  the  steady  plane  flow  problem  described  here  with  the  truncated  usecond-order 
fluid”  constitutive  equation  9,10  w’e  produce  a  solution  very  close  to  that  obtained  when 
the  full  constitutive  equation  is  used  (when  the  flow  is  sufficiently  slow).  By  the  Tanner- 
Giesekus  theorem,  the  Stokes  solution  for  the  velocities,  u,  with  a  modified  pressure,  p, 
satisfies  the  second-order  fluid  equations  [9,10;.  The  resulting  stress  field  is 


a--  \p-  T(uip„  +  u2p,2):  I  r?A 


-  r)T{ B  -  A2  -  [U2I)  -r  t?7”(A2  -  -y2I) 


where  A  and  B  are  the  first  two  Rivlin-Ericksen  tensors;  r),T  and  T'  are  the  constants  of 
the  second-order  fluid  expansion,  p  is  the  Stokes-flow  pressure  field,  and  -y2  =  <rA2 /2. 


j ^  _  *ul,l  ul,2  +  u2,l 

U 1 12  f  “2,1  “2UJJ 


The  symmetries  of  the  Stokes  solution  imply  that  t/i  is  even  and  u2  is  odd  about  <£  as 
functions  of  Zj.  Thus  uji  =  u2>2  =  0  on  <£  (ij  =  0)  and 


A  =  0  'y(xa)' 

[7(12)  0 


Therefore 


A2  -  72I 


and  the  T'-term  in  eq.  (10)  vanishes  on  the  slot  centerline  for  a  second-order  fluid.  The 
r)T- term  is 


-»T(B  -  -1*1) 


B  =  —A  •*?■*  (Vu)^A 


—  U2-4t;.2 


_  u  2Uj.n  (U].2  -  U2.l),, 

1  (**1,2  +  V  2,l)i  j  ~  2t/  j .  1  1 

We  observe  first  that  uii2J  =  «i,i2  =  0,  and  since  u2  is  an  odd  function  about  <£,  u2.n  =  0, 
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We  do  not  believe  that  it  can  be  argued  that  uin  is  zero  on  £.  It  may  be  that  ui,n  is 
“small”  (numerical  evidence  suggests  this),  but  we  do  not  believe  it  can  be  entirely  ignored; 
we  know  that  u^j  =  uliU1  =  0,  and  thus 

Ui(6xj,X2)  =  uj(0, X2)—  -(6xj)-u,.n  ~ 

Were  u^u  zero  or  negligibly  small,  this  would  imply  a  flatness  to  the  streamlines  near  <£ 
which  does  not  seem  to  correspond  to  experiment  or  computation. 

So  we  retain  uj  n  and  observe  that 


Ayu  = 
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and  thus  on  <£ 


a  =-(p-  r«ip„)I-i-  r\ 
-VT 
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From  eq.  (13),  we  may  draw  several  important  conclusions: 

I.  On  £ 

I  o  12  =  vi 

\  °\\  ~  °-n  =  2r?T(72  -  2u1u1>n) 

II. 


1 


) 


P'  =  ~P~  2nT^w^  =  N\JA 


(13) 


(14) 


(15) 


where  estimate  II  applies  when  d  is  such  that  ti,  =  0  or  ulfl ,  =  0  (tij  -  0  at  the  top  of  the 
die  always).  This  will  be  the  case  for  d  =  oc  and  approximately  so  for  d  »  b.  Numerical 
results  suggest  that  d,  b  >  1  is  enough  to  make  Pe  =  Ar)'u,/4  a  good  approximation,  and 
this  in  turn  is  close  to  Niw  4  of  undisturbed  flow  if  h/b  >  2  |4.11  . 


4.  Consequences  of  Tanner-Pipkin  Results 

First  we  observe  that  some  of  the  previous  developments  seem  to  correspond  very 
closely  to  computed  results  in  slow  flows  over  transverse  slots.  Px  is  substantial  and 
positive,  in  spite  of  the  disturbance  induced  by  the  slot,  if  h  b  =  2  and  d  b  =  4.  In  Table  1, 
computations  with  a  Maxwell  fluid  at  De  =  0.25  are  tabulated.  We  find  that  P,  =  N?w/6 
(see  Figure  3).  This  is  slightly  lower  than  Aiu,  6.  but  still  significant.  Evidently  the 
addition  of  the  slot  cannot  have  contributed  much  to  f ^  pd\  .  since  we  have  nearly  the 
same  thrust,  P  =  -0221  at  the  top  of  the  die  as  we  would  have  without  the  slot.  The  only 
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FIGURE  3 

Location  of  undisturbed  and  disturbed  values  (indicated  by 
superscript  zero)  of  quantities  important  to  the  present  work. 


4 


*  *1 


yt 

if 

VI 


'y,‘ 


Nlw/6\N?J*  Pi 


i  0.25  !  1062  I  1008  |  88.33!  80.00  181.12 


TABLE  1 

Disturbed  and  undisturbed  values  of  quantities  in  Figure  3  for 
a  Maxwell  fluid  at  a  low  De  (probably  in  the  second-order  fluid 
range  [7]). 


effect  the  small  disturbance  at  the  opposite  wall  seems  to  have  is  to  lower  the  shear-rate, 
and  thus  .Vj.  at  the  opposite  wall.  The  only  effect  this  has  on  P  there  is  to  lower  it  to 
one  sixth  of  .Y[JU..  as  it  would  be  in  undisturbed  Poiseuille  flow  with  the  slightly  lower 
shear-rate. 

Table  2  shows  the  same  calculations  with  a  modified  Johnson-Segalman  fluid  [12]  at 
De  =  0.24.  ATiu.  and  ou  differ  so  significantly  from  the  Maxwell  case,  because  a  large 
retardation  time.  A,  necessitates  a  higher  shear  rate  to  achieve  a  similar  De  -  (T  -  A)^. 


TABLE  2 

Disturbed  and  undisturbed  values  of  quantities  in  Figure  3  for 
a  modified  Johnson-Segaiman  ;12i  fluid  at  a  low  Dt. 


Here,  too,  we  find  that  P\  = 

For  h/b  >  2,  the  disturbance  due  to  the  slot  has  the  effect  of  deforming  the  streamlines 
in  Stokes  flow  downward  towards  the  slot  mouth  slightly  (see  Figure  4  of  ref.  12,  for 
example).  An  effect  of  this  is  that,  while  012(0,0)  =  0  in  undisturbed  flow,  0i2(O,  i'2)  -  0 
in  disturbed  flow,  for  x2  <  0  w'hich  has  small  magnitude  with  respect  to  h.  Next  we 
observe  that  if  h/b  is  sufficiently  large,  the  disturbance  to  channel  flow  induced  by  the  slot 
will  not  be  of  sufficient  intensity  to  disturb  o 22(0,  x2)  significantly  from  the  value 

^22(0,  h/2)  —  -Pi  =  ~Niw/6 

We  note  that  in  undisturbed  plane  Poiseuille  flow,  022 (^’>*2)  is  constant  in  x 2.  On  the 
other  hand,  by  Tanner  and  Pipkin’s  results  P\  -  P2  =  N®w/4,  and  thus  - P2  =  022(0*  ~d~ 
hj 2)  =  i\\‘w 12.  A  value  near  this  is  confirmed  numerically  for  0  <  De  <  0.3.  Putting 
this  all  together,  we  see  that  022(0,12)  changes  sign  from  negative  to  positive  as  we  move 
from  the  wall  opposing  the  slot  towards  the  bottom  of  the  slot,  and  it  changes  sign  on  the 
slot  side  of  the  channel  centerline,  at  a  point  significantly  closer  to  the  slot  than  (0,  x2).  In 
sum,  observation  of  Stokes  flow  solutions,  the  implications  of  these  observations  in  eq.  (13), 
and  the  Tanner-Pipkin  hole-pressure  result  lead  us  to  deduce  three  important  qualitative 
features  of  the  second-order  fluid  flow  (at  least  for  sufficiently  large  h/b),  illustrated  in 
Figure  4: 

(1)  The  existence  of  x2  <  0  with  small  magnitude  with  respect  to  h. 

(2)  022(O,i2)  is  significantly  negative. 

(3)  022(0.  - d  -  hi 2)  is  significantly  positive. 

5.  HPBL  in  Stress  Integral  Form 

We  refer  to  eq.  (la)  as  the  stress  integral  form  of  the  HPBL  relation  (creeping  flow 
case).  From  the  discussion  above  and  the  derivation  of  ref.  1.  it  is  clear  that  a  more  precise 
form  of  the  relation  should  take  the  disturbance  at  the  wall  opposing  the  slot  into  account. 
This  results  in 

„  r~  a’i  , 


*22  =  -Nfj  6 


(0,ii) 


°22{0,x'2)  <  0 


X„_., 


0-22  —  +^lu >/  12 


FIGURE  4 

Qualitative  picture  of  flow  over  a  transverse  slot  in  second- 
order  fluid  range.  Important  features  are  (a)  that  o 22  is  sig¬ 
nificantly  negative  at  mid-channel  and  in  the  upper  half  of  the 
channel,  but  (b)  o22  does  change  sign  somewhere  below  mid¬ 
channel  and  above  the  bottom  of  the  slot. 


We  should  emphasize  that  our  experience,  guided  by  experimental  and  numerical  evidence 
2.3,6— 8 , ,  is  that  eq.  (16)  is  a  remarkably  accurate  prediction  of  Pe ,  even  in  a  range  of 
flow-rates  far  beyond  the  slow  flows  considered  here.  We  contend,  however,  that  there  are 
several  puzzling  flaws  in  the  derivation  of  eq.  (16),  at  least  when  applied  to  Poiseuiile  flow. 


5.1  Observation  on  HPBL  derivation: 


The  ba^ic  form  derived  by  Higashitani  and  Pritchard  is 


f  V2 

Jy-A 


A’,  dr 


-d-h/2)  2t  dVi 


where  “2/2”  >s  a  streamline  coordinate  (taken  as  a  function  of  the  Cartesian  coordinates 
in  the  integration  limits).  The  assumptions  under  which  the  derivation  is  carried  out  (see 


t 


Appendix  II)  imply  the  limits  may  be  any  stream  function  values  on  (£  and 


o 22(0,X^)  -  <722(0, r£)  = 


dr 


(18) 


where  x \  =  22(1/2)  an<^  x2  ~  £2(1/2)- 


5.2  Three  Flaws  in  HPBL 

Flaw  A:  Piecewise  application  is  contradictory.  Eq.  (18)  is  inconsistent  with  the 
qualitative  picture  of  the  observed  features,  (l)-(3),  illustrated  in  Figure  4,  which  may  be 
further  abstracted  from  any  specific  numerical  results,  as  follows 
Assume: 

(i)  The  existence  of  x'2,  where  <712(0,  x'2)  =  0. 

(ii)  <722(0,  h/2)  has  some  fixed  value  and  <722(0,  h /2)  •  <*22(0,  £•>)  >0. 

(iii)  <722(0,  h/2)  •  022(0,  -d  -  h/2)  <  0. 

Then  rechanging  the  variable  in  eq.  (16) 


Pe 


dy-i 


r~N±d  rV2lk/t)  AT1  dr 
Jo  2r  2r  dt/2 

rvs  (A/2)  y 

/  —  5—  dy2  =  022(0,12)  -  o22 (0,/i/ 2) 

2r  <7j/2 


(19) 


The  crucial  step  follows  from  the  fact  that  the  variable  change  can  make  no  distinction 
between  t/2(~d  -  hi 2)  and  J/2 (^2 ) ’  since  both  locations  o12  =  0.  But  the  qualitative 
picture  abstracted  in  (i) — (iii)  yields 


P(  =  032(0. h/2)  -  022(0,  -d  -  A/2)i  >  1022(0. /t/2)  -  022(0, Xj)|  • 


The  contradiction  lies  in  the  fact  that  the  inequality  is  strict  by  (i) — (iii),  yet  eq.  (19) 
argues  both  sides  should  be  equal.  In  fact,  on  the  basis  of  (1)  —  (3).  one  might  expect  that 
the  inequality  should  be  satisfied  by  a  wide  margin,  since  on  the  left,  the  stresses  add,  but 
on  the  right,  two  nearly  equal  stresses  are  subtracted.  Numerical  evidence  shows  that  the 
two  sides  of  the  above  inequality  are  not  close. 

Flaw  B:  On  —  022  is  not  viscometric  on  (£.  This  seems  to  invalidate  the  change  of 
variable  required  to  derive  eq.  ( 16)  from  eq.  (17)  and  perhaps  more  disturbing,  implies  that 
the  integral  in  eq.  (17)  may  not  even  exist.  The  change  of  variable  requires  that  A  i  can 
be  written  as  a  unique  function  of  r  =  2)-  This  is  not  possible  in  general  as  eq.  (14) 
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shows  that  -  022  is  multiple- valued  as  a  function  of  012  on  <£.  The  multiple-valuedness 
can  be  seen  as  follows: 

Observation:  Consider  that  point  x2  £  (  h  2.h,2)  where  012  -  0.  There  Ui(x2)  *  0,  and 
tii.uf-To)  ^  0,  and  thus  at  x2. 


0; 


0 11  0'1'> 

2t)T 


2u\U 


i“i,n 


f  0 


On  the  other  hand,  for  a  deep  enough  hole.  t>n  0 22  —  0  with  012  ■-  0.  We  thus  have 


0U  -  022 


{0. 


4r]Tu\U  1,11  at  0, 


at  012  -  0; 
also  at  0j2  =  0. 


(20) 


Numerical  results  strongly  support  these  conclusions.  In  fact,  at  (0,x'2)  where  012  =  0, 
and  0n  022  v  0,  there  is  a  pole  in  (01 1  -  022)/'t.  These  poles  have  been  observed 
experimentally  by  I).  C.  Baird  3  .  Note  that  we  have  refrained  from  referring  to  011  -  022 
as  U.Y|,"  reserving  that  notation  for  the  viscometric  function.  What  we  are  saying  here  is 
that,  on  <£,  the  stress-difference  is  not  the  viscometric  function  it  is  presumed  to  be  by  the 
HPBL  derivation,  because  the  flow  is  evidently  not  one  of  constant  curvilinear  shearing 
history.  Thus  01 1  --  022  does  not  take  on  the  values  it  would  have  in  a  viscometric  flow 
with  the  same  A  which  occurs  at  any  given  point  on  the  <£.  Another  way  of  saying  this 
is  that  in  the  <£  flow  A  is  of  viscometric  form,  but  ^  is  not  zero-  as  it  would  be  in  a 
viscometric  flow.  Evidence  suggests  that  uii/i,n  is  quite  small  in  most  flows  where  the 
assumptions  about  A  on  the  are  approximately  valid.  The  troublesome  problem  is  that, 
even  if  u  ] u  1,1 1  is  very  small,  whenever  it  is  not  identically  zero,  but  -j  -  0,  the  centerline 
integral  of  HPBL  does  not  seem  to  exist. 

But  in  Poiseuille  flow  there  appears  to  be  an  even  more  serious  flaw.  Discussion  of 
Flaw  C  is  prefaced  by  the  observation  that  the  assumptions  of  HPBL  should  apply  to 
undisturbed  Poiseuille  flow  (d  =  0).  In  this  case,  ujUj  n  4  0  is  not  a  problem,  and  A’| 
and  0 1 2  have  viscometric  values,  but  HPBL  is  contradictory:  by  eq.(16) 


022(0. /i/2)  —  <^22(0,0) 


■  n: 


di 


■V. 

1 


(21) 


but  this  thrust  difference  should  clearly  be  zero.  Technically.  Higashitani  and  Pritchard  s 
derivation  reduces  to  “0  ~  0V  for  undisturbed  flow  and  needs  a  more  complicated  argument 
to  cover  the  case  in  which  the  streamlines  are  not  curved.  However,  we  believe  that  this 
problem  is  only  technical,  and  the  real  problem  is  that  HPBL  assumes  away  the  streamwise 
gradient,  k. 

Flaw  C:  The  derivation  ignores  the  driving  gradient.  The  equation  in  streamline 
co-ordinates  which  is  integrated  from  point  x2  to  i2  to  give  the  thrust  difference  is  not 


1  dSu  1  dS  12  2  5i2  dS22  _0 

/ll  dj/i  ^  /l 2  dj/2  ^2  5n  -  522  C>y2 

where  Sl:  are  physical  components  of  stress  in  streamline  coordinates.  This  follows  from  the 
geometric  assumptions  of  Higashitani  and  Pritchard,  which  imply  hi  and  h2  =  unknown, 
h3  =  1  and  |^  =  0  unless  t  =  1  and  j  =  2;  and 


[  2  1 _ dhi 

l  11  J  h\dy2' 

(  1  )  j  1  1  =  1  dhi 
1  21  J  \  12  J  hidy2 


The  critical  equation  of  motion  is  then  the  first,  which  has  the  added  term  that 

Higashitani  and  Pritchard  assume  to  be  zero.  In  undisturbed  flow  this  term  is  am,  but  eq. 
(9)  shows  On  i  =  in  undisturbed  plane  Poiseuille  flow.  According  to  eq.  (16)  and  our 
observation  about  eq.  (19),  were  the  HPBL  derivation  correct,  a  substantial  contribution 
to  the  integral  must  come  from  the  upper  half  of  the  channel,  where  the  flow  is  nearly 
undisturbed,  and  the  pressure  gradient  is  nearly  its  undisturbed  value.  It  therefore  seems 
unlikely  that  can  be  ignored.  Leaving  this  term  in  and  continuing  as  Higashitani 

and  Pritchard  do  implies  the  result  of  Flaw  C. 

The  derivation  of  Appendix  II  would  need  some  extra  arguing  to  apply  directly  to 
the  case  of  undisturbed  flow  but  appears  to  be  valid  in  that  case  nevertheless,  since  the 
inclusion  of  the  missing  term  resolves  the  anomoiy: 


1  ft/2(h/2)  ft  j  qt 

<*22(h  12)  -  o22(0)  =  --  /  — (rv- 

2  Ju^to)  T  ”2  vy 2 


k)h2dy2 


And  in  this  flow 


1  dr 

h2dy2  °U'2 


k  =  0 


From  which  the  desired  result  follows.  Note  that  the  key  to  the  validity  of  the  result  is 
including  the  term  omitted  in  the  derivation  of  ref.  1. 

Except  in  the  undisturbed  plane  Poiseuille  case,  the  new  integral  does  not  seem 
amenable  to  a  change  of  variable  to  shear  stress.  Furthermore,  the  presence  of  the  mate¬ 
rial  derivative  terms,  the  poles,  multiple-valued  nature  of  the  integrand  seem  to  render  it 
intractable,  or  even  nonexistent.  The  best  that  can  be  said  is  that  for  a  second  order  fluid, 
the  relation  reduces  to  -N j/4. 


6.  Conclusions 

The  HPBL  relation  does  not  seem  to  be  applicable  to  plane  Poiseuille  flow.  Its  predic¬ 
tion  seems  correct  in  a  second  order  fluid,  but  even  there  its  derivation  is  suspect.  However, 
a  subsequent  report  [7]  will  demonstrate  that  numerical  results  uphold  the  relation 

Ni  =  dPe 

2  a  da 

as  a  reasonably  accurate  measurement  relation.  This  is  the  relation  with  practical  mea¬ 
surement  value,  and  though  it  can  be  formally  derived  via  HPBL,  the  authors  believe  an 
independent  verification  (possibly  numerical)  or  analytic  derivation  would  provide  a  more 
convincing  testimony  to  its  worth  than  the  current  derivation. 
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APPENDIX  I  :  NOTATION 

A  -  First  Rivlin-Ericksen  tensor,  matrix  notation. 

B  -  Second  Rivlin-Ericksen  tensor,  matrix  notation. 
b  -  Slot  width. 

C  -  Constant  of  integration. 

<£  -  Centerline  of  slot. 
d  -  Slot  depth. 

De  -  Deborah  number  =  (T  —  A)^. 

-  Material  derivative. 

/  -  Arbitrary  scalar  function. 
g  -  Arbitrary  scalar  function. 
h  -  Die  height. 

hi  -  Scale  factors  for  streamline  co-ordinates. 

I  -  Identity  tensor,  matrix  notation. 
k  -  Streamwise  pressure  gradient. 

L  -  Length  of  fl  in  ii  -direction. 

A  -  Retardation  time. 

N 1  ~  an  -  ^22  =  Sll  —  ^22- 
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-  Ni  at  wall  in  undisturbed  flow. 

-  Ni  at  wall  opposing  slot. 

-  <722  —  <*33  =  ^22  —  ^33 

-  Hydrostatic  pressure  function. 

-  p  in  Stckes-flow. 

-  Thrust,  -<722- 

-  P  at  wall  opposing  slot. 

-  P  at  slot  bottom. 

-Pi  ~  P2- 

-  L/2  —  6/2. 

-  Stress  tensor  in  streamline  system  (physical  components). 

-  Typical  relaxation  time. 

-  Second-order  fluid  parameter  determining  N? / Ni . 

-  Velocity  vector,  rheological  co-ordinates. 

-  ui  at  channel  center  in  undisturbed  flow. 

-  Rheological  co-ordinates. 

-  Points  on  <£. 

-  Streamline  co-ordinates. 

-  Streamline  co-ordinates  of  x%  and  x|. 

-  Partial  differentiation  with  respect  to  streamline  co-ordinates. 
-±(trA2/2)*. 

-  7  at  wall  in  undisturbed  flow. 

-  7  at  wall  opposing  slot. 

-  Viscosity. 

-  First  normal-stress  coefficient. 

-  Problem  domain,  truncated  channel  of  length  L ,  with  or  without  a  slot. 

-  Stress  tensor,  rheological  co-ordinates, 

-  <7i2  =  S 12  on  <£. 

-  <7 12  at  wall  in  undisturbed  flow. 

-  r  at  wall  opposing  slot. 

-  Christoffel  symbols  for  streamline  co-ordinates. 

-  Partial  differentiation  with  respect  to  rheological  co-ordinates. 
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APPENDIX  II  :  MODIFIED  HPBL  DERIVATION 

Equations  of  motion  in  streamline  coordinates 


ha  d  .hihjka ^ 
hihzhadyj-  hahj  1 


■Sik=0  a  =1,2 


H-P  assumption:  Jj*  =  0  except  when  *  =  l,j  =  2.  By  definition 


ff]  = 


h?  0  0 

0  hi  0 
0  0  1 


{ }  2j.v  ( 


d9<rf)  +  d£  a  <r 

dya  dy0 


d£ai 
dy*  , 


For  a  =  1 


For  a  —  2 


f  2  1  =  dhi 
1 11  /  hi  dy2 


i  ah! 
hi  dy2 


+  {  21  }  h^ki  ^21  ^  {  12  }  h'h2S' 

_  1  3Sn  1  3521  2  dhi 

hi  3yi  h2  dy2  hih2  dy2  21 


_h2_  _3_  hih2  \  ,  /  2  \  h2 
hih2  dyy  h2hyJ  +  \  11  J  h?  11 

i.l  ^  /hi  .  hi  dhi  h2  c 

hi  dyi {5l2)  +  hi  dy2  h2  ^  ^  dy2  h?  5u 


1  dhi  1  dS22 

■ - -S22  H - — 

hih2  dy^  hi  dy2 


1  dhx 
hih2  dy2  11 


and  finally 


Integrating  along  <£ 


1  dhi 
hi  dy2 

*5*1 1  —  $22 
25i2 


dya 

(5n  —  S22) 


_l_dS2i  \ 
h2  #y2  / 


_/*  (±*^  +  2.5*^ 

Jy*  \^i  dyi  h2  dy2  /  25ia 


Note  that  in  undisturbed  flow  on  (£  is  the  Cartesian  covariant  derivative  so  that 
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